Correlations between the nuclear breathing mode energy and properties of 

asymmetric nuclear matter 



(N 

o 

C 

m 



Lie- Wen Chen 1,2 and Jian-Zhong Gu 3 

department of Physics, Shanghai Jiao Tong University, Shanghai 200240, China 
2 Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 
3 China Institute of Atomic Energy, P. 0. Box 275(10), Beijing 102413, China 

(Dated: January 4, 2013) 

Based on microscopic Hartree-Fock + random phase approximation calculations with Skyrme 
interactions, we study the correlations between the nuclear breathing mode energy Sisgmr and 
properties of asymmetric nuclear matter with a recently developed analysis method. Our results 
indicate that the -Eisgmr of 208 Pb exhibits moderate correlations with the density slope L of the 
symmetry energy and the isoscalar nucleon effective mass mJ ? o besides a strong dependence on the 
incompressibility Ko of symmetric nuclear matter. Using the present empirical values of L = 60 ±30 
MeV and mJ,o = (0-8 ± 0.1)m, we obtain a theoretical uncertainty of about ±16 MeV for the 
extraction of Ko from the -Eisgmr of 20S Pb. Furthermore, we find the Sisgmr difference between 
100 Sn and 132 Sn strongly correlates with L and thus provides a potentially useful probe of the 
symmetry energy. 

PACS numbers: 21.65.-f, 24.30.Cz, 21.60.Jz, 21.30.Fe 
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I. INTRODUCTION 

Determination of the equation of state (EOS) for 
isospin asymmetric nuclear matter (ANM) is among fun- 
damental questions in both nuclear physics and astro- 
physics. Knowledge on the nuclear EOS is important 
for understanding not only the structure of finite nu- 
clei, the nuclear reaction dynamics, and the liquid-gas 
phase transition in nuclear matter, but also many crit- 
ical issues such as properties of neutron stars and su- 
pernova explosion mechanism in astrophysics In 
the past more than 30 years, significant progress has 
been made in determining the EOS of symmetric nuclear 
matter from subsaturation density to about 5 times nor- 
mal nuclear matter density pq by studying the nuclear 
isoscalar giant monopole resonances (ISGMR) , collec- 
tive flows Q and subthreshold kaon production [8|, Q in 
nucleus-nucleus collisions. On the other hand, the isospin 
dependent part of the nuclear EOS, characterized essen- 
tially by the nuclear symmetry energy E sym (p), is still 
largely uncertain [H @ . Lack of knowledge on the sym- 
metry energy actually hinders us to extract more accu- 
rately the EOS of symmetric nuclear matter. Therefore, 
to explore and narrow down the uncertainties of both 
the theoretical methods and the experimental data is of 
crucial importance for extracting more stringently infor- 
mation on the nuclear EOS. 

During the past more than 30 years, it has been estab- 
lished that the nuclear ISGMR provides a good tool to 
probe the nuclear EOS around the nuclear normal den- 
sity. In particular, it is generally believed that the in- 
compressibility Kq of symmetric nuclear matter can be 
extracted from a self-consistent microscopic theoretical 
model that successfully reproduces the experimental IS- 
GMR energies as well as the ground state binding en- 
ergies and charge radii of a variety of nuclei [Io| . Ex- 
perimentally, thanks to new and improved experimen- 



tal facilities and techniques, the ISGMR centroid energy 
-EisgmRj i-e., the so-called nuclear breathing mode en- 
ergy, of 208 Pb (a heavy, doubly- magic nucleus with a 
well-developed monopole peak) has been measured with 
a very high precision (less than 2%). Indeed, a value 
of ^isgmr = 14.17 ± 0.28 MeV was extracted from 
the giant monopole resonance in 208 Pb based on an 
improved a-scattering experiment (another value of 
#isgmr = 13.96 ± 0.20 MeV was extracted in Ref. [13). 
The -Eisgmr of 208 Pb has been extensively used to con- 
strain the Kq parameter in the literature 0, [Tl] - |2ll | . It is 
thus important to estimate and eventually narrow down 
the theoretical uncertainty of extracting Kq from the nu- 
clear ISGMR. Theoretically, in fact, it has been found 
that the uncertainty of the density dependence of the 
symmetry energy has significantly influenced the precise 
extraction of the Ko parameter from ISGMR in 208 Pb 
and it also provides an explanation for the observed 
model dependence of the Kq extraction from the ISGMR 
in 208 Pb based on non-relativistic and relativistic mod- 
els 



[14lll5l|22H24|. 

In the present work, we estimate the theoretical un- 
certainty when one extracts the Kq parameter from the 
nuclear ISGMR based on microscopic Hartree-Fock (HF) 
± random phase approximation (RPA) calculations with 
Skyrme interactions. In particular, we study the correla- 
tions between the ISGMR centroid energy and properties 
of ANM with a recently developed analysis method 25| in 
which instead of varying directly the 9 parameters in the 
Skyrme interaction, we express them explicitly in terms 
of 9 macroscopic quantities that are either experimentally 
well constrained or empirically well known. Then, by 
varying individually these macroscopic quantities within 
their known ranges, we can examine more transparently 
the correlation of the ISGMR centroid energy with each 
individual macroscopic quantity and thus estimate the 
theoretical uncertainty of the ISGMR centroid energy 



2 



based on the empirical uncertainties of the macroscopic 
quantities. Our results indicate that the density slope L 
of the symmetry energy and the isoscalar nuclcon effec- 
tive mass m* can significantly change the -Eisgmr. of 
208 Pb and the present uncertainties of L and m* can 
lead to a theoretical uncertainty of about ±16 MeV for 
the extraction of K . Wc further find the -Eisgmr differ- 
ence between 100 Sn and 132 Sn displays a strong correla- 
tion with L and thus provides a potential probe of the 
symmetry energy. 



II. METHODS 

A. Skyrme-Hartree-Fock approach and 
macroscopic properties of asymmetric nuclear 
matter 

The EOS of isospin asymmetric nuclear matter, given 
by its binding energy per nucleon, can be expanded to 
2nd-ordcr in isospin asymmetry 5 as 



spin matrices. The a, £o — £3, xq — X3 are the 9 Skyrme 
interaction parameters which can be expressed analyti- 
cally in terms of 9 macroscopic quantities po, E (p ), the 
incompressibility Kq , the isoscalar effective mass m* , 
the isovector effective mass m* , E sym (po), L, the gradi- 
ent coefficient Gs, and the symmetry-gradient coefficient 
G v [HH3, i- e -> 
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where p = p n + p P is the baryon density with p n and 
p p denoting the neutron and proton densities, respec- 
tively; S = (p rl — Pp)/(p P + p n ) is the isospin asymmetry; 
Eq(p) = E(p,S = 0) is the binding energy per nucleon 
in symmetric nuclear matter, and the nuclear symmetry 
energy is expressed as 



l d*E{p,5) , 
E sym (p) - 1«=0- 
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Around po, the symmetry energy can be characterized by 
using the value of E sym (po) and the density slope param- 
eter L = Zp dEsY a ™ {p) \ p=P0 , i.e., 



L ,p- po, 



E sym (p) = E sym ( Po ) + + Q{{ P -^f). (3) 
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In the standard Skyrme Hartrec-Fock approach, the 
nuclear effective interaction is taken to have a zero-range, 
density- and momentum-dependent form (26| . i.e., 

Via(R,r) = £ (1 + x Pa)S(r) 



£ 1 (l + x 1 P CT )(K' 2 ( 5(r) + ( 5(r)K 2 ) 



+ t 2 (l+x 2 P a )K -6(r)K 

+ iWoi^ + a 2 ) ■ {K x 6(r)K], (4) 

with r = ri— r 2 and R = (ri+r 2 )/2. In the above expres- 
sion, the relative momentum operators K = (Vi— V 2 )/2i 
and K' = — (Vi — V 2 )/2i act on the wave function on the 
right and left, respectively. The quantities P a and Uj de- 
note, respectively, the spin exchange operator and Pauli 



where fc° = (1.57r 2 p ) 1/3 , ^sym(Po) = E syia ( Po ) - 
^sym(Po) — D, and the parameters C, D, a, j3, 7, and 
y are defined as [28| 
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kin 



sym(Po) — 

Bm ( ^"^o ) ■ 1 R tne arj0ve J the isoscalar effective mass 
m* o is the nucleon effective mass in symmetric nuclear 
matter at its saturation density po while the isovector 
effective mass m* corresponds to the proton (neutron) 
effective mass in pure neutron (proton) matter at baryon 
number density p (See, e.g., Refs. [2(1 HU). Fur- 
thermore, Gs and Gv are respectively the gradient and 
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symmetry-gradient coefficients in the interaction part of 
the binding energies for finite nuclei defined as 

£grad = G S {V P f/{2p) - Gy [V(p„ - Pp )] 2 /(2p). (20) 

B. HF + continuum-RPA calculations 

Since the energy of the giant monopolc resonance is 
above the single particle continuum threshold, a proper 
calculation should, in principle, involve a complete treat- 
ment of the particle continuum. In the present work, we 
study the ISGMR of nuclei by using a microscopic HF 
+ continuum-RPA calculations with Skyrme interactions 
[29| . The RPA response function is solved in the coordi- 
nate space with the proton-neutron formalism including 
simultaneously both the isoscalar and the isovector cor- 
relation. In this way, we can take properly into account 
the coupling to the continuum and the effect of neutron 
(proton) excess on the structure of the giant resonances 
in nuclei near the neutron (proton) drip lines [29j . 

The RPA strength distribution of ISGMR of nuclei 

s(e x ) =J2\< n \Q\° > \ 2s ( E * - K «) ( 21 ) 

n 

can be obtained by using the isoscalar monopole operator 
gA=o,-c 1 ^2 (22) 

Furthermore, one can define the k-th energy moment of 
the transition strength by 

m k = J dE x E k x S{E x ). (23) 

The average energy of ISGMR can be defined by the ratio 
between the moments mi and mo, i.e., 

Eave = mi/ mo. (24) 

In addition, the so-called scaling energy of ISGMR can 
be expressed as 

Esca = \Zm 3 /m 1: (25) 

while the ISGMR energy obtained from the constrained 
HF approach can be written as 

E con = \Jm\jm-\. (26) 

The ISGMR energies defined by Eqs. f24 >(p6]l will be- 
come identical in the case of a sharp single peak exhaust- 
ing 100% of the sum rule. In practice, it is found that 
both the experimental data and the theoretical calcu- 
lations show a large width of a few MeV even in the 
most well-established ISGMR in 208 Pb. However, it is 
interesting to note that E ave and E con are rather close 



within a 0.1 ~ 0.2 MeV difference even when the IS- 
GMR peak has a large width although the scaling en- 
ergy E sca has a large uncertainty due to the high en- 
ergy tail of monopole strength, which is always the case 
in experimental data (and see the theoretical results 
in the following). Furthermore, from the relation of 
the energy moments mk+iTrik-i > one can obtain 
E sca > E ave > E con . Therefore, the average energy 
E ave is usually defined as the ISGMR centroid energy 
and compared between the experimental data and the 
theoretical calculations. It should be noted that the sit- 
uation of light nuclei may be quite different from that of 
medium and heavy nuclei considered in the present work 
since the strength distribution of ISGMR for light nuclei 
is usually very fragmented [30I432T ] . It was suggested in a 
recent work [33j that the fragmentation of the strength 
distribution for the light nuclei might be explained by the 
clustering effects which are not considered in the present 
work. 



III. RESULTS 

In the present work, as a reference for the correlation 
analyses performed below with the standard Skyrme in- 
teractions, we use the MSL0 parameter set [25[, which is 
obtained by using the following empirical values for the 
9 macroscopic quantities: po = 0.16 fm~ 3 , Eq( P q) = —16 
MeV, K = 230 MeV, m* = 0.8m, m* v0 = 0.7m, 
E sym {p ) = 30 MeV, L = 60 MeV, Gy = 5 MeV-fm 5 , 
and Gs = 132 MeVfm 5 . And the spin-orbit coupling 
constant Wo = 133.3 MeV -fm 5 is used to fit the neutron 
Pi/2 — P3/2 splitting in 16 O. It has been shown [25| that the 
MSL0 interaction can describe reasonably the binding en- 
ergies and charge rms radii for a number of closed-shell or 
semi-closed-shell nuclei. We further find that the MSL0 
parameter set predicts a value of 1.06 MeV for the split- 
ting of the neutron 3p shell in 208 Pb, which reasonably 
describes the experimental value of 0.9 MeV. It should 
be pointed out that the MSL0 is only used here as a ref- 
erence for the correlation analyses. Using other Skyrme 
interactions obtained from fitting measured binding en- 
ergies and charge rms radii of finite nuclei as in usual 
Skyrme paramctrization will not change our conclusion. 

As numerical examples, in the present work, we choose 
the spherical closed-shell doubly-magic nuclei 208 Pb, 
100 Sn, and 132 Sn. Thus, we do not include the pair- 
ing interaction since it has negligible effects on these 
spherical closed-shell doubly-magic nuclei considered in 
this work (34|. In addition, the two-body spin-orbit and 
the two-body Coulomb interactions are not taken into 
account in the present continuum-RPA calculations al- 
though the HF calculations include both of the inter- 
actions. As pointed out in Ref. (35[, the net effect of 
the two interactions in RPA decreases the centroid en- 
ergy of ISGMR in 208 Pb by about 300 kcV. It should be 
stressed that, in the present work, we do not intend to 
extract accurately the value of the Kq parameter by com- 
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FIG. 1: (Color online) The ISGMR energies of 208 Pb obtained 
from SHF + continuum- RPA calculations with MSL0 by vary- 
ing individually L (a), G v (b), Gs (c), E (po) (d), £ ; sym (p ) 
(e), K (f), m* j0 (g), m* >0 (h), p (i), and W (j). The three 
lines from upper to lower in each panel correspond to E aca , 
Eave, and E con , respectively. 



paring the measured ISGMR ccntroid energy with that 
from HF + continuum-RPA calculations, and our main 
motivation is to explore the theoretical uncertainty for 
extracting K . Meanwhile, we are mainly interested in 
the ISGMR centroid energy difference between 100 Sn and 
132 Sn rather than their respective absolute value. There- 
fore, we do not expect that the two-body spin-orbit and 
the two-body Coulomb interactions in RPA will signifi- 
cantly change our conclusion and further work is needed 
to see how exactly the two interactions in continuum- 
RPA calculations will affect our results. Furthermore, in 
the following calculations, the sum rules are obtained 
by integrating the RPA strength from excitation energy 
E x = 5 MeV to E x = 35 MeV in Eq. 



Isospin scalar giant monopole resonance in 

208p b 



Shown in Fig. Q] arc the ISGMR energies, i.e., 
E sca , E ave , and E con of 208 Pb obtained from SHF + 
continuum-RPA calculations with MSL0 by varying in- 
dividually L, G v , G s , E (p ), E sym (p ), K 0> 
m* , Po, and Wo, namely, varying one quantity at a 
time while keeping all the others at their default val- 
ues in MSL0. Firstly, one can see clearly the ordering of 
E sca > E ave > E con as expected. In particular, for the 
default parameters in MSL0, we obtain E sca = 14.962 
MeV, E ave = 14.453 MeV, and E con = 14.338 MeV. We 
note that the centroid energy of ISGMR E ave = 14.453 
MeV is essentially in agreement with the measured value 
of 14.17 ± 0.28 MeV for the ISGMR in 208 Pb (a 
more recent experimental value of 13.96 ± 0.20 MeV 
was extracted in Ref. Ill)- The agreement will be- 
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come much better if the two-body spin-orbit and the two- 
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FIG. 2: (Color online) SHF + continuum-RPA response func- 
tions of 208 Pb with Skyrme interaction MSL0 by varying in- 
dividually Ko (a), L (b), and m* (c). 



body Coulomb interactions are taken into account in the 
continuum-RPA calculations since the net effect of the 
two interactions in RPA reduces the centroid energy of 
ISGMR in 208 Pb by about 300 keV [Hj]. These features 
imply that the default Skyrme parameter set MSL0 can 
give a good description for the ISGMR in 208 Pb. Further- 
more, one can see from Fig. [T]that, within the uncertain 
ranges considered here for the macroscopic quantities, 
the ISGMR energies display a very strong positive corre- 
lation with Kq as expected. On the other hand, however, 
the ISGMR energies also exhibit moderate negative cor- 
relations with both L and m* while weak dependence on 
the other macroscopic quantities. These results indicate 
that the uncertainties of L and to* may significantly in- 
fluence the extraction of Kq by comparing the theoretical 
value of the ISGMR energies of 208 Pb from SHF + RPA 
calculations with the experimental measurements. 

As for the correlation analysis method in Fig. [1] we 
would like to stress that, when the macroscopic quan- 
tities (except Eq(pq) and po) change individually from 
their base values in MSL0 within the empirical uncer- 
tain ranges considered here, the values of the binding 
energy or the charge rms radii of finite nuclei will vary 
by only about 2% at most (See, e.g., Figs. 4 and 5 in 
Ref. [23|). Therefore, the original agreement of MSL0 
with the experimental data of binding energies or charge 
radii of finite nuclei will essentially still hold with the 
individual change of the macroscopic quantities. In this 
way, changing the macroscopic quantities individually in 
the present correlation analysis approach is equivalent 
to constructing a number of different Skyrme interaction 
sets which can give reasonable description on the ground 
state binding energy and charge rms radii of finite nuclei. 
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The key point and the most important advantage of the -280 
present analysis approach is that in the present correla- 
tion analysis, one knows exactly what is the difference 
among different Skyrme interaction sets constructed by 
using different macroscopic quantities. Furthermore, it 
should be mentioned that the centroid energy of ISGMR 
in heavy nuclei are not sensitive to the values of Eq(pq) 
and po as shown in Fig. [TJ and thus in principle we can 
adjust Eq(pq) and po to give better description for the 
ground state binding energy and charge rms radii of fi- 
nite nuclei without changing significantly the results of 
ISGMR and thus the conclusions in the present work. 

In order to see the dependence of the detailed struc- 
ture of ISGMR in 208 Pb on the values of Kq, L and m* , 
we show in Fig. [2] the SHF + continuum- RPA response 
functions of 208 Pb with MSLO by varying individually 

K , L, and m* „, i.e., Kq = 200 and 270 MeV, L = 30 Qne can expresg K ^ N> z) as 
and 90 MeV, and m* = 0.6m and 0.9m. As can be 




60 90 
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0.7 0.8 
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FIG. 3: The K sa ,t,2 parameter obtained from SHF with MSLO 
by varying individually L (a) and mt,o 3 )- 



L s,0 

seen in Fig. [2l the RPA result displays a single collective 
peak in each case, consistent with the experimental data 
[ 1 lL f36| . Furthermore, it is seen that varying the value of 
K from 200 MeV to 270 MeV strongly shifts the single 
collective peak from about 13.3 MeV to 15.4 MeV while 
varying the value of L (m* ) from 30 MeV (0.6m) to 90 
MeV (0.9m) shifts the single collective peak from about 
14.6 (15.0) MeV to 13.9 (13.9) MeV. These results are 
consistent with the results shown in Fig. [TJ In addition, 
the calculated width with MSLO by varying individually 
Kq, L, and m* shows roughly the same value as that 
of experimental data [Til . l36l |. This agreement implies 
that the width of ISGMR is essentially determined by 
the Landau damping and the coupling to the continuum, 
which are properly taken into account in the present cal- 
culations. 

The ISGMR energy -Eisgmr is conventionally related 
to a finite nucleus incompressibility Ka{N,Z) for a nu- 
cleus with N neutrons and Z protons (A = N + Z) by 
the following definition 



-'ISGMR 



<H 2 Ka(N, Z) 
m (r 2 ) 



(27) 



where m is the nucleon mass and (r 2 ) is the mean square 
mass radius of the nucleus in the ground state. Similarly 
to the semi-empirical mass formula, the finite nucleus 
incompressibility K A (N,Z) is usually expanded as [l(| 

K A (N,Z) = K Q + K smi A- 1/3 +K CUIV A- 2 / 3 



+(K T + K SS A- 1 / 3 ) ( 



N-Z 
A 



+K C 



oul 



(28) 



Neglecting the curvature term K curv , the surface sym- 
metry term K ss and the other higher-order terms of the 
finite nucleus incompressibility Ka{N,Z) in Eq. 



K A (N,Z) = K + K suli A^ 3 + K T 
Z 2 

+K Co 



N-Z 
A 



A 4 / 3 '' 



(29) 



where Kq, K SUT {, K T , and i^ cou i represent the volume, 
surface, symmetry, and Coulomb terms, respectively. 
The K T parameter is usually thought to be equivalent 
to the isospin dependent part, i.e., the K sat ,2 parame- 
ter, of the isobaric incompressibility coefficient of ANM 
(incompressibility evaluated at the saturation density of 
ANM) defined as 



K sat {5) = K + K sata 5 2 + 0{6 4 ). 



(30) 



We would like to point out that the K sat ,2 parameter 
is a theoretically well-defined physical property of ANM 
[H, H3 while the value of the K T parameter may de- 
pend on the details of the truncation scheme in Eq. (|2"8"|) 
[38H42| . Here, we assume i^ S at,2 has similar influences on 
Ka(N,Z) as K T and thus K sat 2 will affect the Eisgmr 
through Eq. ()27|) . and then we can analyze the L and 
m* o dependences of -Eisgmr from the correlations of 
K sa t,2 parameter with L and m* . 

The effects of the density dependence of the symmetry 
energy on the ISGMR centroid energy E ave of 208 Pb has 
been extensively investigated in the literatur e fl4l . [l5l [22| — 
[24j| . It was firstly proposed by Piekarewicz |22| that the 
different symmetry energies used in the non-relativistic 
models and the relativistic models may be responsible for 
the puzzle that the former predicted an incompressibility 
in the range of Kq = 210 — 230 MeV while the latter pre- 
dicted a significantly larger value of Kq w 270 MeV from 
the analysis of the ISGMR centroid energy. It is seen 
from Fig. [JJthat a larger L value (as in usual relativistic 
models) leads to a smaller E ave value and thus a larger 
Kq value is necessary to counteract the deceasing of E ave 
due to a larger L value. Furthermore, Fig. [JJ shows 
that E ave displays a very weak dependence on -E sym (po), 
which is in contrast to the results in Ref. 15( where E ave 
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FIG. 4: The a parameter obtained from SHF with MSLO by 
varying individually Kq (a), Eo(po) (b), m* (c), and po (d). 



is shown to be sensitive to E sym (po). This is due to the 
fact that a constrain on the value of E sym (p = 0.1 fm~ 3 ) 
was imposed in Ref. [l5[, which leads to a strong linear 
correlation between -E S ym(po) and L as shown recently in 
Ref. [H]. 

The symmetry energy dependence of the ISGMR cen- 
troid energy of 208 Pb can be understood from the fact 
that the ISGMR in 208 Pb does not constrain the compres- 
sion modulus of symmetric nuclear matter but rather the 
one of neutron-rich matter, i.e., the isobaric incomprcss- 
ibility coefficient in Eq. (f30|) . From Eq. ([30)) it is clear 
that the ISGMR in 208 Pb (with an isospin asymmetry 
of S = 0.21) should be sensitive to a linear combination 
of Kq and K sat ,2- The if S at,2 parameter is completely 
determined by the slope and curvature of the symmetry 
energy at saturation density as well as the third deriva- 
tive of the EOS of symmetric nuclear matter (see, e.g., 
Ref. [Ill)- Fig- [3] shows the if S at,2 parameter from SHF 
with MSLO by varying individually L and to* . As can 
be seen in Fig. [31 the if sa t,2 parameter decreases with 
both L and m* , and thus K A (N,Z) for 208 Pb will de- 
crease correspondingly if the if sa t,2 parameter has simi- 
lar effects on Ka(N, Z) as the K T parameter and the K ss 
term as well as the other higher-order terms in Eq. ([28]) 
are not important for Ka(N, Z). These results provide 
an explanation on the behavior that the ISGMR energies 
decrease with L and to* observed in Fig. [TJ 

To understand more clearly why the ISGMR energies 
decrease with m* observed in Fig. [TJ it is useful to 
note the fact that, with the standard Skyrme interaction, 
the K and to* cannot be chosen independently if the 
Skyrme interaction parameter a in Eq. (j?]), Eq(pq) and 
po are fixed j44|. It should be stressed here that, instead 



of assuming a fixed value of a as in the usual paramctriza- 
tion and correlation analysis [H],[26j], in the present work, 
the a parameter is determined by four macroscopic quan- 
tities, i.e., Kq, Eq(pq), to* and po as shown in Eq. (fTF 



and thus Kq and m* can be chosen independently. Ne- 
glecting the isospin dependence (assuming N sa Z), the 
nuclear breathing mode energy for medium and heavy 
nuclei can be approximated by j44| 



E 



ISGMR 



'ft 2 (A"o - 63a) 
to (r 2 ) 



{Kq in MeV). 



(31) 



Eq. (|31j) implies that the nuclear breathing mode energy 
can be closely related to both Kq and m* if the parame- 
ter a is free and the values of Eq(pq) and pq are fixed. In 
Fig. El we show the a parameter obtained from SHF with 



MSLO by varying individually Kq, Eq(pq) 



's,0> 



and pq. 



One can see clearly that the a parameter indeed exhibits 
a strong correlation with Kq as expected. However, it 
also displays a moderate dependence on to* , a small 
dependence on Eq(pq), and a very weak correlation with 
Pq. As can be seen in Fig. [H the a parameter increases 
with m* , leading to smaller ISGMR energies according 
to Eq. (j§Tj). which is consistent with the results shown 
in Fig. [JJ In addition, the fact that i^ S at,2 parameter de- 
creases with to* observed in Fig. [SJwill also be partially 
responsible for the behavior of ISGMR energies decreas- 
ing with to* as seen in Fig. [TJsince a smaller if S at,2 value 
will lead to a smaller Eisgmr as discussed previously. 

The above results indicate that the ISGMR centroid 
energy of 208 Pb exhibits moderate correlations with both 



L and 



's,0 



besides a strong dependence on Kq. The ac- 



curate knowledge on L and to* is thus important for 
a precise determination of the Kq parameter from the 
ISGMR centroid energy of 208 Pb. In recent years, sig- 
nificant progress has been made in determining L and 
its value is essentially consistent with L = 60 ± 30 MeV 
depending on the observables and methods used in the 
studies gl [H-i|. Using L = 60 ± 30 MeV, we can 
estimate an uncertainty of about ±0.281 MeV for the IS- 
GMR centroid energy in 208 Pb from Fig. [JJ On the other 
hand, for the isoscalar effective mass, the empirical value 



of; 



(0.8 ± 0.1)to has been obtained from the anal- 



ysis of both isoscalar quadrupole giant resonances data 
in doubly closed-shell nuclei and single-particle spectra 
[44l I56l - l59l ]. From Fig. [1] we can obtain an uncertainty 
of about ±0.382 MeV for the ISGMR centroid energy in 
208 Pb using the empirical value of to* = (0.8 ± 0.1)to. 
Assuming the two uncertainties due to the present uncer- 
tainties of L and to* on the ISGMR centroid energy in 
208 Pb arc independent, we thus can add them quadrati- 
cally to obtain an uncertainty of about ±0.474 MeV for 
the ISGMR centroid energy in 208 Pb. Then, using the 
approximate relation (6K /K ) = 2(<5-Bi S gmr/-Eisgmr) 
from Eq. ([2"T)l . we can obtain an uncertainty of ±7% for 
Kq with -Ejsgmr ~ 14 MeV, namely, about ±16 MeV for 
Kq = 230 MeV. 

Furthermore, including other uncertainties due to Gy , 
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G s , E (p ), E sym (p ), to* , p and W with empiri- 
cal values of G v = ± 40 MeV, G s = 130 ± 10 MeV, 



E {po) 



-16 ± 1 MeV, E sym (p ) = 30 ± 5 MeV, to* 



(0.7 ± 0.1)m, po = 0.16 ± 0.01 fm" 3 and W = 130 ± 20 
MeV, and assuming all the uncertainties are indepen- 
dent, we can obtain from Fig. [1] a total uncertainty of 
about ±0.647 MeV for the ISGMR centroid energy in 
208 Pb, which gives an uncertainty of about ±9% for Kq, 
namely, about ±21 MeV for K Q = 230 MeV. 



B. Isospin scalar giant monopole resonances in 
100 Sn and 132 Sn 



To see the isotopic dependence of the ISGMR centroid 
energy, we study here the spherical closed-shell doubly- 
magic nuclei 100 Sn and 132 Sn. Shown in Fig. [5] are the 
ISGMR centroid energy E ave of 100 Sn and 132 Sn obtained 
from SHF + RPA calculations with MSL0 by varying in- 
dividually L, G v , G s , Eq(po), E sym (p ), K Q , m* , m* , 
Po, and Wq. One can see that the results for neutron- 
rich nucleus 132 Sn are quite similar to those for 208 Pb as 
shown in Fig. [TJ On the other hand, for the symmetric 
nucleus 100 Sn, it is interesting to see that the dependence 
of E ave on the isospin relevant macroscopic quantities, 
namely, L, Gv, E synl (po), to* is very weak. We have 
also checked the case of the stable nucleus 90 Zr for the 
correlation analysis as in Fig. [5j and we find the results 
are very similar to the case of 100 Sn, namely, display- 
ing a much weak correlation with the L parameter while 
a stronger correlation with the Gs parameter compared 
with the case of 208 Pb. This may be understandable from 
the fact that the 90 Zr has a smaller isospin asymmetry, 
i.e., (N — Z)/A = 0.11 compared with 208 Pb where we 
have (iV — Z)/A = 0.21. In addition, the surface coeffi- 
cient Gs may become more important for lighter nuclei 
as expected, leading to a stronger correlation with the 
Gs parameter. From these results, it seems that the 
ISGMR of a heavier and more symmetric nucleus, where 
the symmetry energy effects will be reduced significantly, 
may be more suitable for extracting the Kq parameter. 
In addition, the different E ave -m* correlations between 
100 Sn and 132 Sn observed in Fig. [5] can be understood 
from the fact that K satl 2 parameter decreases with m* 
as shown in Fig. [3j leading additional decrement of E ave 
with to* for the neutron- rich nucleus 132 Sn. 

It is instructive to see the ISGMR centroid energy dif- 
ference between 100 Sn and 132 Sn, which is shown in Fig. 
[6] with MSL0 by varying individually L, Gv, Gs, Eq(pq), 
Esym(po), I<o, m* S fl, m* fi , p , and W . It is very interest- 
ing to see from Fig. [6] that, within the uncertain ranges 
considered here for the macroscopic quantities, the IS- 
GMR centroid energy difference displays a very strong 
correlation with L. However, on the other hand, the IS- 
GMR centroid energy difference exhibits only moderate 
correlations with to* and to* while weak dependence 
on the other macroscopic quantities. These features im- 
ply that the ISGMR centroid energy difference between 
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FIG. 5: (Color online) Same as Fig. [Q but for the ISGMR 
centroid energy E ave of 100 Sn and 132 Sn. The results of 100 Sn 
shift down by 1.5 MeV for a more clear comparison with those 
of 132 Sn. 
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FIG. 6: Same as Fig. [Qbut for the ISGMR centroid energy 
difference between 100 Sn and 132 Sn. 



100 Sn and 132 Sn provides a potential probe of the L pa- 
rameter. Furthermore, it is seen that the ISGMR cen- 
troid energy difference displays opposite correlation with 
to* q and to* q , namely, increases with to* while de- 
creases with to* . Recently, a constraint of to* — to* = 
(0.126 ± 0.051)to has been extracted from global nucleon 
optical potentials constrained by world data on nucleon- 
nucleus and (p, n) charge- exchange reactions [54| . Im- 
posing the constraint m* — to* = (0.126 ± 0.051)to, 
we can expect from Fig. [6] that the correlation of the 
ISGMR centroid energy difference with to* and m* 
will become significantly weak, making the ISGMR cen- 
troid energy difference really a good probe of the L pa- 
rameter. Our results indicate that a precise determina- 
tion of the ISGMR centroid energy difference between 
100 Sn and 132 Sn will be potentially useful to constraint 
accurately the symmetry energy, especially the L pa- 
rameter. This provides strong motivation for measur- 
ing the ISGMR strength in unstable nuclei, which can 
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be investigated at the new/planning rare isotope beam 
facilities at CSR/HIRFL and BRJF-II/CIAE in China, 
RIBF/RIKEN in Japan, SPIRAL2 / GANIL in France, 
FAIR/GSI in Germany, and FRIB/NSCL in USA. 



IV. SUMMARY 

The isoscalar giant monopole resonances of hnite nu- 
clei have been investigated based on microscopic Hartree- 
Fock + random phase approximation calculations with 
Skyrme interactions. In particular, we have studied the 
correlations between the ISGMR centroid energy, i.e., the 
so-called nuclear breathing mode energy, and properties 
of asymmetric nuclear matter within a recently devel- 
oped correlation analysis method. Our results indicate 
that the ISGMR centroid energy of 208 Pb displays a very 
strong correlation with K Q as expected. On the other 
hand, however, the ISGMR centroid energy also exhibits 
moderate correlation with both L and to* while weak 
dependence on the other macroscopic quantities. Using 
the present empirical values of L = 60 ± 30 MeV and 
to* = (0.8 ± 0.1)m, we have obtained an uncertainty 
of about 0.474 MeV for the ISGMR centroid energy 
in 208 Pb, leading to a theoretical uncertainty of about 
±16 MeV for the extraction of Kq from the -Eisgmr. of 
208 Pb. Including additionally other uncertainties due to 
G v , Gs, E (po), E aym (p ), m* , p and W with em- 
pirical values of G v = ± 40' MeV, G s = 130 ± 10 
MeV, Eo(po) = -16 ± 1 MeV, E sym (p ) = 30 ± 5 
MeV, to* = (0.7±0.1)m, p = 0.16 ± 0.01 fin" 3 and 
Wq = 130 ± 20 MeV, we have estimated a total uncer- 
tainty of about ±21 MeV for the extraction of Kq by 



assuming all the uncertainties are independent. These 
results show that the accurate knowledge on L and m* 
is important for a precise determination of the Kq param- 
eter by comparing the measured ISGMR centroid energy 
of 208 Pb with that from Hartree-Fock + random phase 
approximation calculations. 

Furthermore, we have investigated how the ISGMR 
centroid energy difference between 100 Sn and 132 Sn corre- 
lates with properties of asymmetric nuclear matter. We 
have found that the ISGMR centroid energy difference 
between 100 Sn and 132 Sn displays a strong correlation 
with the L parameter while weak dependence on the 
other macroscopic quantities. This feature implies that 
the ISGMR centroid energy difference between 100 Sn and 
132 Sn provides a potentially useful probe of the nuclear 
symmetry energy. Our results also provide strong mo- 
tivation for measuring the ISGMR strength in unstable 
nuclei, which can be investigated at the new/planing rare 
isotope beam facilities around the world. 
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